Suppression of timing errors in short overdamped Josephson junctions 
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The influence of fluctuations and periodical driving on temporal characteristics of short over- 
damped Josephson junction is analyzed. We obtain the standard deviation of the switching time 
in the presence of a dichotomous driving force for arbitrary noise intensity and in the frequency 
range of practical interest. For sinusoidal driving the resonant activation effect has been observed. 
The mean switching time and its standard deviation have a minimum as a function of driving fre- 
quency. As a consequence the optimization of the system for fast operation will simultaneously lead 
to minimization of timing errors. 
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The Rapid Single Flux Quantum (RSFQ) electronic 
devices are promising candidates to built a petaflop com- 
puter due to high operating frequencies of RSFQ ele- 
ments close to 1 THz 0. Moreover they are of particu- 
lar interest in solid-state quantum information process- 
ing. They may be used both for realization of qubits and 
for characterization of a macroscopic quantum behavior, 
e.g., readout electronics for quantum computing 2]. The 
processes going on in RSFQ devices are based on a re- 
production of quantum pulses due to spasmodic changing 
by 2ir of the phase difference of damped Josephson junc- 
tions (JJ). The major restriction in the development of 
RSFQ logic circuits is given by influence of fluctuations 
III liil Idl LJ ■ Moreover, the operating temperatures of the 
high-T c superconductors lead to higher noise levels by 
increasing the probability of associated errors. Recently 
lots of investigations were performed to study pulse jitter 
and timing errors in RSFQ circuits 0, . The Timing 
errors is one of reasons limiting a 16-bit RSFQ micropro- 
cessor prototype clock frequencies to 20GHz instead 
of theoretically predicted hundreds GHz. Therefore, in- 
vestigation of possible ways for suppression of time jitter 
of transmitting signals is a very important problem from 
fundamental and technological points of view. An RSFQ 
circuit in fact has long timing loops, so the time jitter 
produces a noise-induced digital accumulative effect in 
Josephson transmission lines 30- It is also important 
for different branches of physics, where a nonlinear ele- 
ment is driven either from an external source or from an- 
other element, as in neural networks. For RSFQ circuits 
three different types of digital errors may be identified. 
First, storage errors may occur during the passive stor- 
age of data: fluctuations can induce a 2-7r phase flip and 
thus switch the quantizing loop into the neighboring flux 
state. The analytical description of the mean time of such 
noise-induced flips, valid for arbitrary noise intensity, has 
been presented in jy]. Second, decision errors may be 
produced in the two junction comparator. This situation 
has been studied in detail in the linearized overdamped 
J J model • If an RSFQ circuit operates at high speed, 
another type of noise-induced error becomes important. 
Due to noise the time interval between input and output 



pulses fluctuates, producing "timing" error. Timing er- 
rors have also been studied in the linearized overdamped 
J J model |J| . It is known that in noisy nonlinear systems 
with a metastable state, the resonant activation (RA) 
HHE3 and noise enhanced stability (NES) [E E3 phe- 
nomena may be observed. These effects may play positive 
and negative role in accumulation of fiuctuational errors 
in RSFQ logic devices: the RA phenomenon minimizes 
timing errors, while the NES phenomenon increases the 
switching time. These effects however were not still ob- 
served in previous investigations due to the use of lin- 
earized models 0, . Moreover the limiting frequencies 
of RFSQ devices and possible optimizations, in order to 
increase working frequencies and reduce timing errors in 
RSFQ circuits, is an open question. 

This letter is aimed to answer this question by in- 
vestigating nonlinear noise properties of an overdamped 
JJ, subjected to periodic driving, to understand possible 
ways of RSFQ circuits optimization for high-frequency 
operation with minimal timing errors. To this end we 
consider the dynamics of a short overdamped J J, under 
a current /, given by the Langevin equation Q 



_-ydip(t) du(ip) 
dt dip 



i F (t), 



(1) 



u(cp) = 1 - cos</>- i(t)(p, i(t) =i + /(*), (2) 

where (p is the order parameter phase difference, u((p) 
is the dimensionless potential profile, f(t) is the driv- 
ing signal, i = I/I c with I c the critical current, and 
ip (t) ~ Ip/I c . Ip is the thermal noise current with 
(I F (t)) = 0,and(I F (t)I F (t + T)) = 2j5(t)/uj c . Here 
7 = 2ekT / (Til c ) is the dimensionless noise intensity, 
T is the temperature, k is the Boltzmann constant, 
lu c = 2eR]\[I c /h is the characteristic frequency of the JJ, 
Rn is the normal state resistance of a J J, e is the elec- 
tron charge and h is the Planck constant. Let, initially, 
the JJ is biased by a current smaller than the critical 
one (io < 1), and the junction is in the superconductive 
state. The current pulse f(t), such that i(t) > 1, switches 
the junction into the resistive state. However, an output 



pulse will be born not immediately, but later on. Such a 
time is the switching time, which is a random quantity 
characterized by mean value and standard deviation. As 
an example of a driving with sharp fronts we consider the 
dichotomous signal f(t) — Asign(sin(arf)), and as an ex- 
ample of a driving with smooth fronts - sinusoidal signal 
f(t) — Asm(u>t). In spite that we consider the peri- 
odic driving, that makes the results more evident, it is 
obvious, that below the cut-off frequency the switching 
occurs during half of the first period, so other periods 
are not important. Besides, our adiabatic analysis may 
easily be generalized to an arbitrary form of f(t). In 
computer simulations we set lu c — 1 and, therefore, in 
plots u is normalized to lu c . The mean switching time 
(MST) t = (t) = J °° tw(t)dt, and its standard deviation 
(SD) <t, may be introduced as characteristic time scales of 
the evolution of the probability P(t) — W(ip, t)dip to 
find the phase within one period of the potential profile 
Eq.© 



w(t) 



dP(tp ,t) 



dt[P(<po,oo)-P(tp ,Q)] 
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We choose therefore tp2 — it, fx = —it and the initial 
distribution at the bottom of a potential well. Let us 



t(co), cr(cJ) 
6.0- 



4.0- 



2.0- 



0.0- 



tM y=0.2 
T (ai), y=0.02 




ajoj), y=0.2 



ajoj), 7=0.02 



0.01 



0.1 



1 

OJ 



FIG. 1: The MST t{lj) and SD a(uj) as functions of the 
frequency for dichotomous driving, for two values of noise in- 
tensity: 7 = 0.2,0.02, and io = 0.5, i = 1.5. The results 
of computer simulations are: t(u>) (solid line) and cr(w) (di- 
amonds and circles). Dashed line are the theoretical results 
given by Eqs.JIjJl, ||H). 

first consider the case of dichotomous driving. From the 
Fokker-Planck equation corresponding to the Langevin 
Eq. JU, we calculate the MST and its SD. In Fig. 1 we 
report the behaviors of MST t(ui) and SD <j(oj) as a func- 
tion of the driving frequency for two values of noise inten- 
sity: 7 = 0.2,0.02. For dichotomous driving both MST 
and SD do not depend on the driving frequency below a 
certain cut-off frequency, above which the characteristics 
degrade. In the frequency range -j- 0.2w c , therefore, we 
can describe the effect of dichotomous driving by time 



characteristics in a constant potential. The exact ana- 
lytical expression of MST and its asymptotic expansion 
are: for arbitrary 7 Q 

e-vMhfy. / e u{ifi)h dif\, (4) 

and for 7 <C 1 

/ x h (<Pa) ~ AOo) +71/2(^2) + /2O0)] + •• 

Tc{V0) = , (5) 

OJc 



2 fitssx(x/2) - 1\ 

arctan | K J ' — | , (6) 



f 2 (x) = l/(2( J -sinx) 2 ). 

Following ref. [lTJ, the exact expression for t 2c — (t 2 ) in 
a time-constant potential may be derived 



T2c(>o) = T c (ip ) - 

2 
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where H(x) = J e u ^'~< J e - u( ^^ J e^^dzdydv, and 

xv y 

T c (ipo) is given by fllf. The asymptotic expression of a = 
y/ Tic — t 2 in the small noise limit 7 <C 1 is 



*(<Po) = — ^27 [%) + /3W] + .., (8) 
F((fo) = /i(V2)/a(^2) - 2fi(<Po)M<Po)+ 
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The comparison between the asymptotic theoretical re- 
sults (Eqs.||SJ,||S}) and simulations is reported in Fig. 1. 
The agreement is very good within the frequency range 
0-^0. 2uj c . It is interesting to see that the SD of the switch- 
ing time scales as square root of noise intensity. The de- 
pendencies of the MST and its SD on the bias current 
are presented in Fig. 2 for 7 = 0.001. The agreement be- 
tween theoretical results and simulations is very good for 
all the bias current values investigated. In the low noise 
limit Eq.© gives actually the same results of ref. 0], 
since the largest contribution to the MST comes from 
the deterministic term T c (ip ) = t\ {h^) — h{Vo)} ■ 
However the Eq. (JSJ, in some cases, significantly devi- 
ates from the linearized calculations. For 7 = 0.001 and 
i = 1.2, a = OAu^ 1 in ref. 4], while we get a = 0.436W" 1 . 
For larger current i — 1.5, the discrepancy is larger: 
a = 0.06^,7 1 in |4|, and we get a — 0. 14w~ . In the 
inset of Fig. 2 we report the behaviour of SD (Eq.©) 
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as a function of noise intensity 7. For noise intensity 
values up to 7 = 0.05 the agreement with computer sim- 
ulations is very good. Therefore not only low temper- 
ature (7 < 0.001), but also high temperature devices 
may be described by Eqs.JHJ), ©■ If noise intensity is 
rather large, the phenomenon of NES may be observed: 
the MST increases with the noise intensity, as it may be 
easily seen from Eq.JSJ. In the design of large arrays of 
RSFQ elements, operating at high frequencies, it is very 
important to consider this effect, otherwise it may lead 
to malfunctions due to the accumulation of digital errors. 
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FIG. 2: The MST and SD vs bias current for the time- 
constant case for 7 = 0.001. Solid lines - formulas @ and 
®, diamonds and circles - results of computer simulation. 
Inset: the SD vs noise intensity for the time-constant case 
(i > 1). Solid line - formula @, circles - results of computer 
simulation for i = 1.5. Dashed line - formula (JHJ, diamonds - 
results of computer simulation for i = 1.2. 

Now let us consider the case of sinusoidal driving. The 
corresponding time characteristics may be derived using 
the modified adiabatic approximation |9lll0| 



P(ip ,t) = exp j-^ 



1 



-dt' 



. „, - , • (9) 

with T c ((po, t') given by l@J). For dichotomous driving the 
value of initial current io has a weak effect on tempo- 
ral characteristics, as in ref. Q, while it is important 
for sinusoidal driving , since it also defines the poten- 
tial barrier height. We focus now on the current value 
i = 1.5, because i — 1.2 is too small for high frequency 
applications. In Fig. 3 the MST as a function of the 
driving frequency for different values of bias current is 
shown. For smaller io the switching time is larger, since 
ip$ = arcsin(io) depends on io- On the other hand, the 
bias current io must be not too large, since it will lead, in 
absence of driving, to the reduction of the mean life time 
of superconductive state, i.e. to increasing storage errors 
(Eq.(0J). Therefore, there must be an optimal value of 
bias current io, giving minimal switching time and ac- 
ceptably small storage errors. Storage errors, in fact, are 




FIG. 3: The MST vs frequency for f(t) = ^sin(wt) (com- 
puter simulations) for i = 1.5. Long-dashed line - 7 = 0.02, 
short-dashed line - 7 = 0.05, solid line - 7 = 0.5, from top 
to bottom io = 0.5; A = 1, io = 0.8; A = 0.7. Inset: com- 
parison between simulations and theoretical results obtained 
from Eq.@ for io = 0.5; A = 1, and 7 = 0.02 (diamonds), 
7 = 0.05 (circles) and 7 = 0.5 (crosses). 



acceptably small up to i = 0.99 for 7 = 0.001 0. We 
observe the phenomenon of resonant activation: MST 
has a minimum as a function of driving frequency. The 
approximation @ does not describe the resonant activa- 
tion effect at high frequencies, but it works rather well 
below 0.1 uj c , that is enough for practical applications. 
Moreover, it is interesting to see that near the minimum 
the MST has a very weak dependence on the noise in- 
tensity, i. e. in this signal frequency range the noise 
is effectively suppressed. We observe also the NES phe- 
nomenon. There is a frequency range, around 0.2 — 0.4oj c 
for io = 0.5 and around 0.3 — 0.5w c for io = 0.8, where 
the switching time increases with the noise intensity. The 
NES effect increases for smaller io because the potential 
barrier disappears for a short time interval within the 
driving period T — 2it/uj |ll| and the potential is more 
fiat, so noise has more chances to prevent the phase to 
move down and delay switching process. This effect may 
be avoided, if the operating frequency does not exceed 
0.2w c . Besides (see Fig. 4) the SD also increases above 
0.2oj c . The plots of SD as a function of driving frequency 
for 7 = 0.02, i = 1.5 and different values of io are shown 
in Fig. 4. The approximation © is not so good for SD 
as for MST, even if the qualitative behaviour of SD is re- 
covered. We see that the minimum of a(oj), for 7 = 0.02, 
is located near the corresponding minimum for t(oj) in 
Fig. 3. For the SD the optimal frequency range, where 
the noise induced error will be minimal, is from 0.1 to 
0.3 for the considered range of parameters. It is interest- 
ing to see that, near the minimum, the SD for sinusoidal 
driving actually coincides with SD for dichotomous driv- 
ing (Eq.(JHJ). This means that even in the case of smooth 
driving, the limiting value of SD may nearly be reached, 
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FIG. 4: The SD vs frequency for /(t) = Asm(ut) and 7 = 
0.02. Computer simulations - dash-dotted line: io = 0.3; A = 
1.2, short-dashed line: io = 0.5; A = 1, long-dashed line: 
io = 0.8; A = 0.7. The MST is given by crosses for comparison 
(io = 0.8; A = 0.7). Formula (g) - solid line. Formula Q - 
diamonds. 



but the RSFQ circuit must be properly optimized. 

Finally we note that close location of minima of MST 
and its SD means that optimization of RSFQ circuit for 
fast operation will simultaneously lead to minimization 
of timing errors in the circuit, which is the main result 
of this paper. 

In the present paper we reported an analytical and 
numerical analysis of influence of fluctuations and peri- 
odic driving on temporal characteristics of the J J. For di- 
chotomous driving the analytical expression of standard 
deviation of switching time works in practically inter- 
esting frequency range and for arbitrary noise intensity. 
For sinusoidal driving the resonant activation effect has 
been observed in the considered system: mean switching 
time has a minimum as a function of driving frequency. 
Near this minimum the standard deviation of switching 
time takes also a minimum value. The RA phenomenon 
was observed very recently in the underdamped Joseph- 
son tunnel junction |l2j . Our theoretical investigation 
could motivate experimental work in overdamped JJs 
as well. Utilization of this effect in fact allows to sup- 
press time jitter in practical RSFQ devices and, there- 
fore, allows to significantly increase working frequencies 
of RSFQ circuits. Our study is not only important to 



understand the physics of fluctuations in a Joscphson 
junction, to improve the performance of complex digital 
systems, but also in nonequilibrium statistical mechan- 
ics of dissipative systems, where noise assisted switching 
between metastable states takes place 0, 0, 0, 0, ^| . 
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